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^2 Abstract. Jeans's (1902 & 1929) gravitational instability criterion gives truly spectac- 
ular errors in its predictions of cosmological structure formation according to Gibson's 
(1996) new theory. It is suggested that the linear perturbation stability analysis leading to 
the acoustic Jeans length scale criterion L > Lj = Vg/ (pG)^^^ is quite irrelevant to grav- 
I itational condensation, which is intrinsically nonlinear and nonacoustic. Instead, conden- 

. . . . 

I sation is limited by viscous or turbulent forces, or by diffusivity, at Schwarz length scales 

o 

Lsv, LsT, or LsD, respectively, whichever is larger, independent of Lj and the sound 

^ \ speed Vs- By these new criteria, cosmological structure formation begins in the plasma 

J> I epoch soon after matter dominates energy, at i « Lsv = il'^/ pG)^^^ scales correspond- 
•i— I , 

ing to protosuperclusters, decreasing to protogalaxies at the plasma-gas transition, where 
7 is the rate-of-strain of the expanding universe, v is the kinematic viscosity, p is the den- 
sity, and G is Newton's gravitational constant. Condensation of the primordial gas occurs 
at mass scales a trillion times less than the Jeans mass to form a 'fog' of micro-brown- 
dwarf (MBD) particles that persist as the galactic baryonic dark matter, as reported by 
Schild (1996) from quasar-microlensing studies. Nonbaryonic dark matter condensation 
is prevented by its enormous diffusivity at scales smaller than Lsd = {D^ / pGY^^ , where 
D is the diffusivity, so it forms outer halos of galaxies, cluster-halos of galaxy clusters, 
and supercluster-halos. 



Key words: Cosmology: theory - dark matter - early Universe - gravitational lensing ■ 
large-scale structure formation of the Universe 



Send offprint requests to: C. H. Gibson 



2 



Carl H. Gibson: The fluid mechanics of dark matter formation 



1. Introduction 

Jeans's theory fails because it is linear. Linear theories typically fail when applied to 
nonlinear processes; for example, applications of the linearized Navier-Stokes equations 
to problems of fluid mechanics give laminar solutions, but observations show that actual 
flows are turbulent when the Reynolds number is large. 

Cosmology (e.g.; Peebles 1993 , Padmanabhan 199^ relies exclusively on Jeans's the- 



ory in modelling the formation of structure by gravitational instability. Because Lj for 
baryonic (ordinary) matter in the hot plasma epoch following the Big Bang is larger than 
the Hubble scale of causality Lh = ct, where c is the speed of light and t is the time, no 
baryonic structures can form until the cooling plasma forms neutral gas. Star and galaxy 
formation models invented to accommodate both Jeans's theory and the observations of 
early structure formation have employed a variety of innovative maneuvers and concepts. 
Nonbaryonic "cold dark matter" was invented to permit nonbaryonic condensations in 
the plasma epoch with gravitational potential wells that could guide the early formation 
of baryonic galaxy masses. Fragmentation theories were proposed to produce MQ-stars 
rather than Jeans-superstars at the IO^Mq proto-globular-cluster Jeans mass of the hot 
primordial gas. 

Both of these concepts have severe fluid mechanical difficulties according to the 



Gibson 1996 theory. Cold dark matter cannot condense at the galactic scales needed 
because its nonbaryonic, virtually coUisionless, nature requires it to have an enor- 
mous diffusivity with supergalactic Lsd length scales. Fragmentation theories (e.g.. 



Low and Lynden-Bcll 19761 ) are based on a faulty condensation premise that implies 



large velocities and a powerful turbulence regime that would produce a first gener- 
ation of large stars with minimum mass determined by the turbulent Schwarz scale 
LsT = e^^^ /{pG)^^^ , where e is the viscous dissipation rate of the turbulence, with a 
flurry of starbursts, supernovas, and metal production that is not observed in globu- 
lar star clusters. The population of small, long-lived, metal-free, globular cluster stars 
observed is strong evidence of a quiet, weakly turbulent formation regime. 



2. Jeans's acoustic theory 

Jeans considered the problem of gravitational condensation in a large body of nearly con- 
stant density, nearly motionless gas. Viscosity and diffusivity were ignored. The density 
and momentum conservation equations were linearized by dropping second order terms 
after substituting mean plus fluctuating values for the density, pressure, gravitational po- 
tential, and velocity. Details of the derivation are given in many cosmological texts (e.g.; 
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Kolb and Turner 1994, p342) so they need not be repeated here. The mean gravitational 
force V(/) is assumed to be zero, violating the Poisson equation 

= 47rGp, (1) 

where (f) is the gravitational potential, in what is known as the Jeans swindle. Cross- 
differentiating the linearized perturbation equations produces a single, second order dif- 
ferential equation satisfied by Fourier modes propagating at the speed of sound Vs . From 
the dispersion equation 

^2 ^ v^k^ - AnGp, (2) 

where u is the frequency and k is the wavenumber, a critical wavenumber kj — 
{4:TrGp/VgY^'^ exists, called the Jeans wavenumber. For k less than fcj, uj is imaginary 
and the mode grows exponentially with time. For k larger than kj, the mode is a propa- 
gating sound wave. Density was assumed to be a function only of pressure (the barotropic 
assumption) . 

Either the barotropic assumption or the linearization of the momentum and density 
equations are sufficient to reduce the problem to one of acoustics. Physically, sound 
waves provide density nuclei at wavecrests that can trigger gravitational condensation if 
their time of propagation X/Vg for wavelength A is longer than the gravitational free fall 
time Tg = (pG)^^/^. Setting the two times equal gives the Jeans gravitational instability 
criterion: gravitational condensation occurs only for A > Lj. 

Jeans's analysis fails to account for the effects of gravity, diffusivity, or fluid me- 
chanical forces upon nonacoustic density maxima and density minima; that is, points 
surrounded on all sides by either lower or higher density. These move approximately 



with the fluid velocity, not Vs, (Gibson 1968). The evolution of such zero gradient 
points and associated minimal gradient surfaces is critical to turbulent mixing the- 
ory ( [Gibson et al. 1988 ). Turbulence scrambles passive scalar fields such as tempera- 
ture, chemical species concentration and density to produce nonacoustic extrema, sad- 
dle points, doublets, saddle lines and minimal gradient surfaces. A quasi-equilibrium 
develops between convection and diffusion at such zero gradient points and minimal 
gradient surfaces that is the basis of a universal similarity theory of turbulent mixing 
( pibson 1991 ) analogous to the universal similarity theory of Kolmogorov for turbulence. 
Just as turbulent velocity fields are damped by viscosity at the Kolmogorov length scale 
Lk = {i^/lY^^, where v is the kinematic viscosity and 7 is the rate-of-strain, scalar fields 
like temperature are damped by diffusivity at the Batchelor length scale Lb = (-D/7)^^^, 
where D is the molecular diffusivity. This prediction has been confirmed by laboratory ex- 
periments and numerical simulations ( Gibson et al. 198^ for the range 10^^ < Pr < 10^, 
where the Prandtl number Pr = v/D. 
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On cosmological length scales, density fields scrambled by turbulence are not nec- 
essarily dynamically passive but may respond to gravitational forces. In the density 
conservation equation 

dp/dt + Vi{dp/dxi) ~ Dcffd^p/dxjdxj (3) 

the effective diffusivity of density Dcs = D — L'^ /tq is affected by gravitation in the 
vicinity of minimal density gradient features, and reverses its sign to negative if the fea- 



ture size L is larger than the diffusive Schwarz scale Lsd (Gibson and Schild 1998a). 
LsD = {P"^ I pG)^l^ is derived by setting the diffusive velocity V£, « D jL of an isodcnsity 
surface a distance L from a minimal gradient configuration equal to the gravitational 
velocity Vg « L/Tg. Thus, nonacoustic density maxima in a quiescent, otherwise ho- 
mogeneous, fluid are absolutely unstable to gravitational condensation, and nonacoustic 
density minima are absolutely unstable to void formation, on scales larger than Lsd- 



Jeans believed from his analysis ( Jeans 1929 ) that sound waves with A > Lj would 



grow in amplitude indefinitely, producing unlimited kinetic energy from his gravitational 
instability. This is clearly incorrect, since any wavecrest that collects a finite quantity of 
mass from the ambient fluid will also collect its zero momentum and become a nonacoustic 
density nucleus. From the enormous Jeans mass values indicated at high temperature, 
he believed he had proved his speculation that the cores of galaxies consisted of hot 
gas (emerging from other Universes!) and not stars, which could only form in the cooler 
(smaller Lj) spiral arms, thrown into cold outer space by centrifugal forces of the spinning 
core. The concepts of pressure support and thermal support often used to justify Jeans's 
theory are good examples of bad dimensional analysis, lacking any proper physical basis. 



3. Fluid mechanical theory 

Gravitational condensation on a nonacoustic density maximum is limited by either dif- 
fusion or by viscous, magnetic or turbulent forces at diffusive, viscous, magnetic, or 
turbulent Schwarz scales Lsx, whichever is largest, where X is D,V,M,T, respectively 
( [Gibson 1996 , Gibson and Schild 1998a ). Magnetic forces are assumed to be unimpor- 
tant for the cosmological conditions of interest. Gravitational forces Fg « p^GL'^ equal 
viscous forces Fy ~ pv^L'^ at Lsv = {^1 1 pGY^'^ , and turbulent forces Ft ~ p{e)'^^^L^^^ 
at LsT = £i/V(pG)3/4. Kolmo gorov's theory is used to estimate the turbulent forces as 
a function of length scale L. 

The criterion (Gibson 1996, 1997a, 1997b; Gibson and Schild 1998a, 1998b) for grav- 
itational condensation or void formation at scale L is therefore 



L>{Lsx)^,^.X^D,V,M,T. 



(4) 
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4. Structures in the plasma epoch 

Without the Jeans constraint, structure formation begins in the early stages of the hot 
plasma epoch after the Big Bang when decreasing viscous forces first permit gravitational 
decelerations and sufficient time has elapsed for the information about density varia- 
tions to propagate; that is, the decreasing viscous Schwarz scale Lsv becomes smaller 
than the increasing Hubble scale Lh = ct, where c is the velocity of light. Low levels 
of temperature fluctuations of the primordial gas indicated by the COsmic microwave 
Background Experiment (COBE) satellite {ST/T « 10^^) constrain the velocity fluc- 
tuations 5v/c ^ 10^^ to levels of very weak turbulence. Setting the observed mass of 
superclusters « 10'*^ kg equal to the Hubble mass pL^ computed from Einstein's equa- 
tions ( Weinberg 1972| , Table 15.4) indicates the time of flrst structure was ~ 10^^ s, or 



30 000 y dGibson 1997b| ). 



Setting Lh « 3 lO^" m (10 kpc) = Lsv gives « 6 10^"^ m^ s^^ with p « lO^i^ kgm^^ 
and J 1/t — 10^^^ rads^^. Such a large viscosity suggests a neutrino-electron- proton 
coupling mechanism, presumably through the Mikheyev-Smirnov-Wolfenstein (MSW) 



effect (Bahcall 1997), supporting the Ncutrino-98 claim that neutrinos have mass. 

The viscous condensation mass pLgy decreases to about 10''^ kg (Gibson 1996) as 



the Universe expands and cools to the plasma-gas transition at t « lO^'^ s, or 300 000 y, 
based on Einstein's equations to determine T and p and assuming the usual dependence of 



viscosity v on temperature T (Weinberg 1972). Assuming gravitational decelerations that 



are possible always occur, we see that protosupercluster, protocluster, and protogalaxy 
structure formation should be well underway before the emergence of the primordial gas. 

5. Primordial fog formation 

The first condensation scales of the primordial gas mixture of hydrogen and helium are 
the maximum size Schwarz scale, and an initial length scale Ljc = {RT/ pG)^^^ equal to 
the Jeans scale Lj but independent of Jeans's linear perturbation stability analysis, and 
acoustics, where R is the gas constant of the mixture. From the ideal gas law p/p — RT 
we see that density increases can be compensated by pressure increases with no change 
in temperature in a uniform temperature gas, and that gravitational forces Fg « p^GL'^ 
will dominate the resulting pressure gradient forces Fp k, pL^ = pRTL^ for length scales 
L > Lie- Taking i? « 5 000 m^s^^K-i, p w 10"!^ kgm"^ ^ ( [Weinberg 1972| ) and 
T Ri 3 000 K gives a condensation mass pL^^j « 10^ Mq, the mass of a globular cluster of 
stars. Because the temperature of the primordial gas was observed to be quite uniform 
by COBE, we can expect the protogalaxy masses of primordial gas emerging from the 
plasma epoch to immediately fragment into proto-globular-cluster (PGC) gas objects on 
Lie ~ 3 10^'' m (10 pc) scales, with subfragments at (isx)niax- 
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The kinematic viscosity v of the primordial gas mixture decreased by a factor of 
about a triUion from plasma values at transition, to i/ « 2.4 10^^ m^ assuming the 
density within the PGC objects are about 10^^^ kgm^'^. Therefore, the viscous Schwarz 
scale Lsv ~ (2.4 lO^^ iq-is/IO^i^ 6.7 lO"")^/^ = 1.910" m, so the viscous Schwarz 
mass Msv ~ L\yp = 6.8 10^^ ^j. ^/^^ ^ g.S lO^"* kg using p = 10~". The turbulent 
Schwarz mass Mst ~ 8.8 10^^ kg assuming 10% of the COBE temperature fluctuations 
are due to turbulent red shifts ([((5u/c)/((5r/r)] 10^^) as a best estimate. 

We see that the entire universe of primordial H-He gas turned to fog soon after the 
plasma-gas transition, with primordial fog particle (PFP) mass values in the range 10^'^ 
to 10^^ kg depending on the estimated density and turbulence levels of the gas. The 
time required to form a PFP is set by the time required for void regions to grow from 
minimum density points and maximum density saddle points to surround and isolate 



the condensing PFP objects (Gibson and Schild 1998a). Voids grow as rarefaction waves 



with a limiting maximum velocity Ys set by the second law of thermodynamics, so the 
minimum PFP formation time is Tppp < {Lsx)ina.x/Vs, or about 10'^ y. Full condensation 
of the PFP to form a dense core near hydrostatic equilibrium requires a much longer time, 
near the gravitational free fall time « 2 10^ y. 

Radiation heat transport during the PFP condensation period before the creation of 
dust should have permitted cooling to temperatures near those of the expanding universe. 
After about a billion years hydrogen dew point and freezing point temperatures (20-13 K) 
would be reached, forming the niicro-brown-dwarf conditions expected for these widely 
separated (10'^ — 10"* AU) small planetary objects that comprise most of the baryonic dark 
matter of the present universe and the materials of construction for the stars and heavy 
elements. Because such frozen objects occupy an angle of less than a micro-arcsecond 
viewed from their average separation distance, they are invisible to most observations 
except by gravitational microlensing, or if a nearby hot star brings these volatile comets 
out of cold storage. 

6. Observations 

A variety of observations confirm the new theory that fluid mechanical forces and diffu- 



sion limit gravitational condensation (Gibson 1996), and confute Jeans's (1902 & 1929) 
acoustic criterion: 



— quasar microlensing at micro-brown-dwarf frequencies (Schild 1996), 

— tomography of dense galaxy clusters indicating diffuse (nonbaryonic) superhalo dark 
matter at Lsd scales with Dni, « lO^^m^s"^ (Tyson and Fischer 1995), 



the Gunn-Peterson missing gas sequestered as PFPs, 
the dissipation of 'gas clouds' in the Ly — a forest. 
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— extreme scattering events, cometary globules, FLIERS, ansae, Herbig-Haro 'chunks', 
etc. 



Evidence that the dark matter of galaxies is dominated by small planetary mass 
objects has been accumulating from reports of many observers that the multiple images 
of lensed quasars twinkle at corresponding high frequencies. After several years spent 
resolving a controversy about the time delay between images Q0957+561A,B, to permit 
correction for any intrinsic fluctuations of the hght intensity of the source by subtraction 
of the properly phased images, Bchild 1996| announced that the lensing galaxy mass 
comprises w IO^^Mq "rogue planets" that are "likely to be the missing mass." 

Star-microlensing studies from the Large Magellanic Cloud have failed to detect lens- 
ing at small planetary mass frequencies, thus excluding this quasar-microlensing popula- 



tion as the Galaxy halo missing mass ( Alcock et al. 1998 , Renault et al. 1998 ). However, 
the exclusion is based on the unlikely assumption that the number density of such small 
objects is uniform. The population must have mostly primordial gas composition since 
no cosmological model predicts this much baryonic mass of any other material, and must 
be primordial since it constitutes the material of construction, and an important stage, 
in the condensation of the gas to form stars. Gravitational aggregation is a nonlinear, 
self-similar, cascade process likely to produce an extremely intermittent lognormal spatial 
distribution of the PFP number density, with mode value orders of magnitude smaller 
than the mean. Since star-microlensing from a small solid angle produces a small num- 
ber of independent samples, the observations estimate the mode rather than the mean, 
resolving the observational conflict ( [Gibson and Schild 1998b|) . 



7. Conclusions 

1. Jeans's gravitational instability criterion L > Lj is irrelevant to gravitational struc- 
ture formation in cosmology and astrophysics, and is egregiously misleading in all of 
its applications. 

2. The correct criterion for gravitational structure formation is that L must be larger 
than the largest Schwarz scale; that is, L > (isx)niax' where X is D, V, M, T, de- 
pending on whether diffusion or viscous, magnetic or turbulent forces limit the grav- 
itational effects. 

3. Structure formation began in the plasma epoch with protosupercluster to protogalaxy 
decelerations. 

4. Gravitational condensations began soon after the plasma-gas transition, forming 
micro-brown-dwarfs, clustered in PGCs, that persist as the dominant dark matter 
component of inner galactic halos (50 kpc). 
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5. The present fluid mechanical theory and its cosmological consequences regarding the 



forms of baryonic and nonbaryonic dark matter (Gibson 1996) is well supported by 



observations, especially the quasar-microlensing of Schild 1996 and his inference that 
the lens galaxy mass of Q0957+561A,B is dominated by small rogue planets (inter- 
preted here as PFPs). 
6. Star-microlcnsing studies that rule out MBDs as the Galaxy missing mass 



( Alcock et al. 1998 , Renault et al. 1998 ), contrary to the quasar-microlensing evi- 
dence and the present theory, are subject to extreme undersampling errors from 
their unwarranted assumption of a uniform number density distribution, rather 
than extremely intermittent lognormal distributions expected from nonlinear ag- 
gregational cascades of such small objects as they form nested clusters and stars 
( [Gibson and Schild 1998b| ). 

Acknowledgements. Numerous helpful suggestions were provided by Rudy Schild. 



References 

Alcock, C, et al. 1998, ApJ, 499, L9 

Bahcall, J. N. 1997, in Unsolved Problems in Astrophysics, Eds. Bahcall, J. N. and Ostriker, J. 

P., Princeton Univ. Press, NJ 
Gibson, C. H. 1968, Phys. Fluids, 11, 2305 

Gibson, C. H., W. T. Ashurst, and A. R. Kerstein 1988, J. Fluid Mech. 194, 261 
Gibson, C. H. 1991, Proc. R. Soc. A, 433, 149 
Gibson, C. H. 1996, Appl. Mech. Rev., 49, 299 

Gibson, C. H. 1997a, The Identification of Dark Matter, Ed.; N. J. C. Spooner, World Scientific, 
114 

Gibson, C. H. 1997b, Dark Matter in Astro- and Particle Physics, Ed.: H. V. Klapdor- 

Kleingrothaus and Y. Ramachers, World Scientific, 409 
Gibson, C. H. and Schild, R. E. 1998a, to be published 
Gibson, C. H. and Schild, R. E. 1998b, to be published 
Jeans, J. H. 1902, Phil. Trans. R. Soc. Lond. A, 199, 1 

Jeans, J. H. 1929, Astronomy and Cosmology, Cambridge Univ. Press, Cambridge, UK 

Kolb, E. W. and M. S. Turner 1994, The Early Universe, Addison- Wesley Publishing Company 

Low, C. and Lynden-Bell, D. 1976, MNRAS, 176, 367 

Padmanabhan, T. 1993, Structure Formation in the Universe, Cambridge University Press, 
Cambridge, UK 

Peebles, P. J. E. 1993, Principles of Physical Cosmology, Princeton Univ. Press, Princeton, NJ 

Renault, C. et al. 1998, A&A, 329, 522 

Schild, R. 1996, ApJ, 464, 125 

Tyson, J. A. and P. Fischer 1995, ApJ, 446, L55 

Weinberg, S. 1972, Gravitation and Cosmology: Principles and Applications of the General 
Theory of Relativity, John Wiley and Sons, New York 



